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Starting from state-by-state calculations of exclusive rates of the ordinary muon capture (OMC), 
we evaluated total /^“-capture rates for a set of light- and medium-weight nuclear isotopes. We 
employed a version of the proton-neutron quasi-particle random phase approximation (pn-QRPA, 
for short) which uses as realistic nuclear forces the Bonn C-D one boson exchange potential. Special 
attention was paid on the percentage contribution to the total -capture rate of specific low-spin 
multipolarities resulting by summing over the corresponding multipole transitions. The nuclear 
method used offers the possibility of estimating seperately the individual contributions to the total 
and partial rates of the polar-vector and axial-vector components of the weak interaction Hamilto¬ 
nian for each accessible final state of the daughter nucleus. One of our main goals is to provide a 
reliable description of the charge changing transitions matrix elements entering the description of 
other similar semileptonic nuclear processes like the charged-current neutrino-nucleus reactions, the 
electron capture on nuclei, the single /3^-decay mode, etc., which play important role in currently 
interesting laboratory and astrophysical applications like the neutrino-detection through lepton- 
nucleus interaction probes, and neutrino-nucleosynthesis. Such results can be also be useful in 
various ongoing muon-capture experiments at PSI, Fermilab, JPARC and RCNP. 
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I. INTRODUCTION 

In recent years, various sensitive experiments take ad¬ 
vantage of the powerful muon beams produced in well 
known muon factories (PSI, Fermilab, JPARG, RGNP in 
Osaka and others) for standard and non-standard muon 
physics probes mm- Among the standard model probes 
those involving muon capture on nuclei specifically those 
emitting X-rays and/or several particles (p, n, a, etc.) 
after /i“-capture, which are important to understand the 
rates and spectra of these particles, are investigated [T]. 
For example, at PSI researchers are interested in exper¬ 
iments based on the emission of charged particles from 
muonic atoms of Al, Si and Ti or neutron emission follow¬ 
ing muon capture from Fe, Ga ,Si and Al [I]. Also very 
recently, in the highly intense muon facilities MuSIC at 
RGNP, Osaka, Japan, nuclear muon capture reactions 
(on Mo, Pb, etc.) are planned in order to study, nuclear 
weak responses (for neutrino reactions, etc.) [ 2 ]. For ex¬ 
periments like the above, it is important, before going to 
the rates of the emitted particles, to know the first stage 
muon capture process. 

As it is well known, when negative muons, yT, pro- 
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duced in a meson factory, slow down in matter, it is pos¬ 
sible for them to be captured in atomic orbits. After¬ 
wards, fast electromagnetic cascades bring these muons 
down to the innermost (Is or 2p) quantum orbits (in this 
way muonic atoms are produced) [3H7]. A bound muon 
in the muonic atom may disappear either by decay known 
as muon decay in orbit or by capture by the nucleus the 
main channel of which is the ordinary muon capture rep¬ 
resented by the reaction [SHin] 

Tb + (^> ^ ~ ^)* + ( 1 ) 

where {A, Z) denotes the initial atomic nucleus with mass 
number A and proton number Z while (A, Z — 1)* stands 
for an excited state of the daughter nuclear isotope. 

The reaction ([^ is a well known example of symbio¬ 
sis of atomic, nuclear and particle physics. In this work, 
however, we will concentrate on its nuclear physics as¬ 
pects. As muon-capture in nuclei presents many advan¬ 
tages for the study of both nuclear structure and the fun¬ 
damental electro-weak interactions UMl] , process 0 
has been the subject of extensive experimental and the¬ 
oretical investigations started early on 50’s using closure 
approximation or sum-over partial rates to find the to¬ 
tal /i“-capture rate (the measured quantity). [3HS1 ITCF 
[ 22 ]. In the plethora of the relevant papers, the most 
important motivation was rested on the hope to explain 
how nucleons (hadronic current) inside the nucleus cou¬ 
ple weakly to the lepton held (leptonic current). The nu- 
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clear physics aspects of process Q , however, still possess 
some yet unresolved fundamental problems, e.g. those 
related to the nucleon-nucleon and lepton-nucleus inter¬ 
actions, the question whether the individual properties 
of the nucleons change when they are packed together 
in the nucleus or remain essentially unaffected like the 
coupling of the nucleon to the leptonic field, etc. 

The interest of studying -capture has recently been 
revived P51 - E5] due to its prominent role in testing the 
nuclear models employed in several physical applications 
in neutrino physics and astrophysics [2HU2H] . Specifically, 
/^“-capture is a very useful test for various nuclear meth¬ 
ods used to describe semi-leptonic weak charged-current 
reactions [sms] as the electron capture in stars (criti¬ 
cal in the collapse of supernovae) [26l [28l [29] , the neu¬ 
trino nucleus scattering (important in the detection of as- 
trophysical neutrinos) [571 HB] , and other reactions m- 
This is due to the fact that the muon-capture involves 
a large momentum transfer and, hence, it can provide 
valuable information about effects which are not found 
in processes like the beta-decay modes (on medium mo¬ 
mentum transfer processes however, useful information 
can be also obtained from low-spin forbidden transitions 
of beta-decays and charge exchange reactions) [3D|. Fur¬ 
thermore, there is an intimate relation between the inclu¬ 
sive muon capture rate and the cross section for the an¬ 
tineutrino induced charged-current reactions, since both 
are governed by the same nuclear matrix elements and 
proceed from the same set of initial to the same final nu¬ 
clear states [nn |5B|-|2H] • Moreover, from the ground state 
transition matrix elements of the /i“-capture process one 
may also derive cross sections for the beta decay modes 
[75] , Calculations on single beta decay which are more 
difficult to calculate, need explicit nuclear structure cal¬ 
culations [ST] , 

The purpose of the present work is to perform detailed 
state-by-state calculations [32H4T| of exclusive muon cap¬ 
ture rates and concentrate on the individual contribu¬ 
tion of each basic multipole operator inducing low-lying 
excitations in the daughter nucleus. In contrast, most 
of the previous muon capture calculations have been 
performed within the assumptions of closure approxi¬ 
mation nzi El mj. Towards this aim, the pn-QRPA 
provides a reliable description of the required nuclear 
transition matrix elements d d [TBHTSl [i5H5D] . Our 
extensive channel-by-channel calculations would be car¬ 
ried out for the exclusive, partial and total muon cap¬ 
ture rates, and the results refer to the nuclear iso¬ 
topes Fe,^^ Zn and which cover 

the light- and medium-weight region of the periodic ta¬ 
ble. We also specialize on the individual contributions of 
the polar-vector and axial-vector components of the /re¬ 
capture operators in each of the multipole states and in 
the total Ordinary Muon Carture (OMC) rate. Despite 
the fact that the semileptonic process Q is studied for a 
long time [Hzl EH EH E9H22I, essentially only the total 
muon-capture rates have been measured for a great num¬ 
ber of nuclear isotopes [ElEnillllS]- On the theoretical 


side, various nuclear methods using several residual in¬ 
teractions allowed the calculation of total capture rates 
on many nuclei with an accuracy of about 10% compared 
to the experimental rates. However, for only few isotopes 
exclusive capture rates to specific states in the daughter 
nucleus have been determined [SlEnilMlES]. As the ex¬ 
perimental data for muon capture rates are quite precise, 
and the theoretical techniques of evaluating the nuclear 
response in the relevant nuclear systems are well devel¬ 
oped 01311 I35j . it is worthwhile to see to what extent 
the exclusive capture rates are theoretically understood. 


Furthermore, we mention that, there appear recently, 
clear indications that the axial-vector coupling constant 
gA = 1.262 in a nuclear medium is reduced from its free 
nucleon value PEniHSlEIHMI- The evidences for such a 
renormalization of the value gA come primarily from the 
analysis of beta decay modes between low-lying states of 
medium-heavy nuclei jSlj but the use of a quenched gA 
value is mainly invoked from the second-order core polar¬ 
ization caused by the tensor force m and the screening 
of the Gamow-Teller (GT) operator by the A-hole pairs 
[56) . Thus, it is necessary to scrutinize on the in-medium 
quenching of the axial vector coupling constant which is 
in agreement with various well-known indications that gA 
is reduced to the value of gA ~ 1.000. In this work we 
are not going to study systematically this effect, but we 
will compare our results of /i-capture rates obtained with 
the values (i) gA = 1.262 and (ii) gA = 1.135 with other 
theoretical ones obtained with the latter value. 


The rest of the paper, is organized as follows. In 
Section 2, we summarize briefly, the main character¬ 
istics of the effective charged-current weak interaction 
Hamiltonian and present the main formalism of the or¬ 
dinary muon capture rates which is based on our com¬ 
pact formalism for the relevant nuclear transition ma¬ 
trix elements (relying on the Donnelly-Walecka projec¬ 
tion method) and in the expressions for exclusive, partial 
and total muon capture rates 0EIII31- Special focus 
is given on the calculation of the nuclear wave functions 
derived within the context of the pn-QRPA. In Section 
3, we concentrate on the determination of the required 
model parameters for the nuclear ground state, derived 
by solving the BCS (Bardeen Cooper Schrieffer) equa¬ 
tions, as well as of the excited states (solution of the 
pn-QRPA equations). Our results (Sestion 4) refer to 
exclusive, partial and total muon capture rates, of the 
above mentioned nuclear isotopes, which cover the light- 
and medium-weight region of the periodic table. We also 
include the individual contributions of the polar-vector 
and axial-vector operators in each of the multipole states 
and in the total Ordinary Muon Carture (OMC) rate. 
Finally, in Section 5, we summarize the main conclusions 
extracted from the present work. 
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II. FORMALISM OF MUON CAPTURE RATES 

The ordinary muon-capture process, that takes place 
in muonic atoms and is represented by the semi- 
leptonic reaction Q, proceeds via a charged-current 
weak-interaction Hamiltonian which is written as a prod¬ 
uct of a leptonic, and a hadronic current, as 

[SmSllMlIST] 

( 2 ) 

where G = GpcosOc with Gp and 9c being the well known 
weak interaction coupling constant and the Cabbibo an¬ 
gle, respectively. 

From the nuclear theory point of view, the main task is 
to calculate the partial and total capture rates of the re¬ 
action 0 which are based on the evaluation of exclusive 
nuclear transition matrix elements of the form 

(/IKN) = J (3) 

(the integration is performed in the region of the nuclear 
system). In the latter expression |i) and |/) denote the 
initial (ground) and the final nuclear states, respectively. 
The quantity stands for the leptonic matrix el¬ 

ement written in coordinate space with q being the 3- 
momentum transfer. The magnitude of is defined from 
the kinematics of the process and is approximately given 
by m 


q = qf = rrifj, - eb + Ei - Ef (4) 


where is the muon rest mass, eb is the muon-binding 
energy in the muonic atom, Ei denotes the energy of the 
initial state of the parent nucleus and Ef the final energy 
of the corresponding daughter nucleus. 

In the unihed description of all semi-leptonic electro- 
weak processes in nuclei developed by Donnelly and 
Walecka nisi [SHE], the calculation of the required tran¬ 
sition strengths of Eq. ^ is based on a multipole decom¬ 
position of the hadronic current density which leads to a 
set of eight independent irreducible tensor multipole op¬ 
erators (four of them come from the polar-vector compo¬ 
nent and the other four from the axial-vector component 
of the nuclear current). In the present work we assume 
that the pn-QRPA excitations \J^) have good quantum 
numbers of angular momentum (J), parity (tt) and energy 
which is a basic assumption for the Donnelly-Walecka 
projection method to be applicable. In this spirit, the 
computation of each partial transition rate of the muon 
capture is written in terms of the eight different nuclear 
matrix elements (between the initial \Ji) and the final 
\Jf) states) as 


Ai->/ — 


2G2g2 


- Cj)U)f (5) 


where represents the muon wave function in the 
Is muonic orbit [44]. The operators in Eq. ([^ refer 
to as Coulomb Mj, longitudinal £j, transverse elec¬ 
tric Tj'' and transverse magnetic multipole op¬ 

erators and contain polar-vector and axial-vector parts 
(see Appendix 0. The factor Rf in Eq. ([^ takes 
into consideration the nuclear recoil which is written as 

Rf = ^I + qf/Mtarg^ , with Mtarg being the mass of 
the target nucleus. 


III. DESCRIPTION OE THE NUCLEAR 
METHOD 

Eor reliable predictions of partial muon-capture rates, 
a consistent description of the structure of the ground 
state \Ji) of the parent nucleus as well as of the multipole 
excitations \Jf) of the daughter nucleus are required. In 
the present work, the state-by-state muon capture rates 
are evaluated using Eq. ([^ with the transition matrix el¬ 
ements between the states \Ji) and \Jf) determined with 
the use of the BCS and pn-QRPA equations, respectively 
(the BCS equations determine the ground state and the 
pn-QRPA equations provide the excited states as it is 
shown below) [34H4T| . To this end, at first we have cho¬ 
sen the active model space (the same for proton and neu¬ 
tron configurations) for each studied isotope consisted of 
the single particle j-shells shown in Table |l| 


TABLE I: The used active model space with the respective 
harmonic oscilator parameter for all the studied nuclei. In 
the last column the major harmonic oscillator shells N plus 
the individual orbits used for each nucleus are listed. 


Nucleus 

b(h.o) 

Model Space 

Core 

Active 

Levels 

N (huj) 


1.809 

No 

10 

0,1,2,3 

32^ 

1.843 

No 

12 

0,l,2,3,0gr9/2,0(/7/2 

48 Ti 

1.952 

No 

12 

0,1,2,3,039/2,0(77/2 


1.996 

16q 

12 

2,3,4 


2.043 

16q 

12 

2,3,4 


2.138 

16q 

16 

2 , 3 , 4 , 0 / 111 / 2 , 0 / 19 / 2 , 1 / 7 / 2 , 1 / 5/2 


As it is well known, in a rather good approximation, 
the nucleus can be considered as a system of Z protons 
and N neutrons moving independently inside the nuclear 
volume and attracted by the nuclear center through a 
central strong nuclear force. This central attraction is 
well described by a mean held which, in our case, is as¬ 
sumed to be a Woods-Saxon potential with a Coulomb 
correction and a spin-orbit parts |3Sj- For the latter po¬ 
tential we tested two different parametrizations: i) that 
of Bohr and Motelson [E], and ii) that of the IOWA 
group |59j and found that both give rather similar re¬ 
sults. For the purposes of the present work, however, we 
adopted the more realistic IOWA parametrization |59j . 

For a reliable nuclear Hamiltonian, in addition to the 
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mean field, the two-nucleon correlations, known as resid¬ 
ual two-body interaction, are necessary to be included. 
Towards this aim, we employed the pn-Bonn C-D one- 
boson exchange potential, but, since the initially eval¬ 
uated bare nucleon-nucleon matrix-elements of the lat¬ 
ter potential refer to all nuclides with mass number A, 
for a specific isotope (A, Z) studied, a renormalization 
of these two-body matrix elements was carried out with 
the use of four multiplicative parameters: The first two, 
known as pairing parameters for protons (p) and 

neutrons (n), renormalize the monopole (pairing) inter¬ 
action which is the part of the correlations involved at the 
BCS level for the description of the considered indepen¬ 
dent quasi-particles. The third, gpp, tunes the particle- 
particle channel and the fourth, gph, renormalizes the 
particle-hole interaction of the Bonn C-D potential. We 
briefly summarize the adjustement of these parameters 
below (subsection III BI. 


A. Determination of the parent nucleus ground 
state 

The ground state of the parent nucleus, is obtained 
within the context of the BCS theory where the one- 
quasi-particle states are deduced by solving (iteratively) 
the BCS equations. Towards this aim one is defining 
quasi-particle creation, a\ and annihilation, a, opera¬ 
tors related to the particle-creation, cj,, and particle an¬ 
nihilation, Ck, operators through the Bogolyubov-Valatin 
transformations [60l |6l] 

al = Ukcj, - UfcC„, = UkC^ + Vkcl, (6) 

where £« denotes the time reversed particle annihila¬ 
tion operator defined as = (—1 )A+™'=c_k with —k = 
(kj—nik)- The probability amplitudes Vk and Uk for the 
k single particle level to be occupied or unoccupied, re¬ 
spectively, are [60] 

2(p,n) _ 1 
- 2 

= 1 — v^) where is the single particle energy of 
the jK-level and Ap (A„) denotes the chemical potential 
for protons (neutrons). Moreover, the solution of the rel¬ 
evant BCS equations gives the single quasi-particle ener¬ 
gies |4S1150] 

£?"> = + (8) 

with being the theoretical energy gaps (A^ = 

~J2k'>o^kkk'k''^k’Vk') [SSj- From the Solution of the gap 
equation [^ ifH] 

Aw = (9) 


1 - 


-A, 


■p{n) 


E- 


,p(n) 


(7) 


(here the notation is, [j] = -^/SJ^PT) one obtains the 
pairing gaps for protons A* and neutrons A^ through 
the renormalization of the proton and neutron pairing 
matrix elements {{kk)0\Q\{k'k')0) of the residual interac¬ 
tion, using the parameters Sp^j^^ g^^ir- The lowest 
quasi-particle energy, obtained from the gap equation, is 
determined, through the pairing parameters gp'^l enter¬ 
ing the theoretical gaps of Eq. ([^ so as to reproduce the 
experimental (empirical) energy gaps Ap^P given from 
the three point formula |46j 


A exp _ 
Pin) - 


1 

4 


5p(„) [(A - 1, Z - 1(Z))] - 25p(„) [(A, Z)] 


+ >5'p(n) [(A-f 1, Z-|-1(Z))] 


( 10 ) 


In the latter equation Sp and Sn are the experimen¬ 
tal separation energies for protons and neutrons, respec¬ 
tively, of the target nucleus (A,Z) and of the neighboring 
nuclei (A ± 1, Z ± 1) and (A ± 1, Z). Here, we used the 
method of Ref. [IS] to obtain the g^’^ir values for the 
studied nuclei and tabulate them in Table El We note 
that, in order to achieve the reproducibility of the exper¬ 
imental energy spectrum in similar QRPA calculations 
some authors modify slightly the Woods-Saxon proton 
and neutron single particle energies in the vicinity of the 
nuclear Fermi surfaces [261136j . In this work, we pay spe¬ 
cial attention on the reproducibility of the energy spec¬ 
trum of the daughter nucleus as is discussed in detail in 
the next section. 


TABLE II: Parameters for the renormalization of the inter¬ 
action of proton pairs, and neutron pairs, Qpair- They 

have been fixed in such a way that the corresponding experi¬ 
mental gaps, Ap^^ and , are quite accurately reproduced. 


Nucleus 

9pair 

jpatr 

^exp ^theor ^exp ^theor 

(MeV) (MeV) (MeV) (MeV) 


1.1312 

1.0601 

3.1428 

3.1429 

3.0375 

3.0377 

325 

0.8862 

0.8230 

2.0978 

2.0979 

2.0387 

2.0386 

48 y. 

0.9259 

0.9833 

1.5576 

1.5578 

1.9112 

1.9111 


0.9866 

0,9756 

1.3626 

1.3626 

1.5682 

1.5683 


1.0059 

0.9271 

1.7715 

1.7716 

1.2815 

1.2814 


0.9057 0.7838 

1.8567 

1.8568 

1.1184 

1.1183 


B. The pn-QRPA excitation spectrum of the 
daughter nucleus 

For the purposes of the present study, transitions be¬ 
tween the |0+) ground state of a rather spherical even- 
even parent-nucleus and the excited states of the re¬ 
sulting daughter nucleus are the basic ingredients. For 
several charged-current reactions, the pn-QRPA method 
provides a reliable description of the nuclear excited 
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states of the resulting odd-odd nuclear system in Eq. Q 
[46j . Here, we exploit this advantage in order to derive 
the excitation spectrum of the daughter nucleus produced 
in the /r-capture process. In this context, we first define 
the two quasi-fermion operators and A (which obey 
boson commutation relations in a correlated RPA ground 
state) as [T^ IBOl IM] 

AUjM) = ( 11 ) 

rrimi'ITT'i) 

ArmiJM) = {-ly-^ArmiJ-M). (12) 

Afterwards, we write down the pn-QRPA phonon opera¬ 
tors 

Q%m = E &^L(^M) + y^,i™( JM)], (13) 

m<i 

V enumerates the multipole states of the multipolarity 
,V, that creates the excitation \v) = | J(f) by acting on 
the QRPA vacuum |0)Qfl;p^ as [T^ 15^ 155115^I55VHT] 

\Ju) = Q'j^mP)QRPA- (14) 

The X (forward) and Y (backward) scattering ampli¬ 
tudes entering Eq. are obtained by solving the pn- 
QRPA equations (pn-QRPA eigenvalue problem) which 
in matrix form is written as |60| 

-;.(?:). (-) 

denotes the excitation energy of the QRPA state 
\jy. Thus, the X and Y amplitudes are calculated seper- 
ately for each multipole set of states (multipolarity). 

The reliability of the QRPA excitations flj„ and of 
the corresponding many-body nuclear wave functions is 
checked through the reproducibility of the energy spec¬ 
trum of the final odd-odd nucleus. The values of particle- 
particle (gpp) and particle-hole {gph) parameters in the 
set of isotopes chosen (determined separately for each 
multipolarity) [351 [351 SI] lie in the region 0.65 - 1.20 
(with the exception of the 1'*" and 2“ multipolarities in 
some isotopes, for which the values are rather small, 0.2 
- 0.6) [62]. Such small values of the strength param¬ 
eters come out in studies of charged current reactions 
(e“-capture, single- and double-beta decays) when fit¬ 
ting simultaneously the QRPA parameter, gpp, and the 
axial vector coupling constant, gA |47| EH IH3H65] . We 
stress that in our QRPA method the strength parame¬ 
ters are determined through the reproduction of the en¬ 
ergy spectrum of the daughter nucleus but we have also 
made an effort to test them through the GT energy po¬ 
sition and the total GT strength [^|35]. Even though 
our GT-type operator contributes differently (due to the 


presence of the Bessel function), we found that, the to¬ 
tal GT strength differs significantly (more than a factor 
of 2.5) from the experimental one, although, the energy 
position is well reproduced. In our muon capture (and 
e“-capture) rates the simultaneous variation of gA and 
gpp parameters has not been checked extensively (see Ref. 

EH- 

We furthermore note that, in order to achieve the re¬ 
producibility of the experimental energy spectrum of the 
daughter nucleus and for measuring the excitation en¬ 
ergies of the daughter nucleus from the ground state of 
the initial (even-even) nucleus, some authors shift the 
entire set of QRPA spectrum by about Xp — A„ in the 
muon capture process [24]. In our present study we also 
adopt the latter treatment, so, the calculated pn-QRPA 
energy spectrum of each individual multipolarity J'^ is 
shifted in such a way that the first calculated value of 
each multipole state (i.e. 1)^, 2(''...etc), to approach as 
close as possible to the corresponding lowest experimen¬ 
tal energy of the daughter nucleus. Such a shifting is 
necessary whenever in the pn-QRPA a BCS ground state 
is used, a treatment adopted by other groups too SZIIIE]. 
Table [ml shows the shifting applied to the QRPA spec¬ 
trum for each multipolarity of the studied nuclei. We 
note that, a similar treatment is required in QRPA cal¬ 
culations for double-beta decay studies where the excita¬ 
tions derived for the intermediate odd-odd nucleus (in¬ 
termediate states) through p-n or n-p reactions from the 
neighboring nuclei do not match each other [47] . The re¬ 
sulting low-energy spectrum (up to 3.0 MeV) using our 
pn-QRPA method, agrees well with the experimental one 
as can be seen from Pig. 

Before proceeding to our results, it is worthwhile to 
briefly summarize the advantages of the calculational 
procedure followed in performing the present detailed cal¬ 
culations of partial and total muon capture rates as com¬ 
pared to the methods used by other groups [8l [T6l [23l l24] . 
In the earlier pioneering work of Foldy and Walecka US], 
the authors related the dipole capture rate to the exper¬ 
imental photoabsorption cross section and used symme¬ 
try arguments to compare polar-vector and axial-vector 
matrix elements. The afore mentioned authors derived 
/i“-capture sum rules based on the GDR strength ex¬ 
cited after /i“-capture. The required GDR amplitudes 
are obtained (for light and medium nuclei) from the cor¬ 
responding photo-absorption cross sections. Later, on 
the calculations of Eramzhyan et al. [23] employed a 
truncated model space with ground state correlations and 
adopted the standard free nucleon coupling constants. In 
the work of Kolbe et al. |5] , for the calculations of muon 
capture rates, use of the continuum RPA method was 
made with the free nucleon form factors, while recently, 
Zinner et al. [23] proposed the use of a quenched value for 
the axial-vector coupling constant gA in order to reliably 
evaluate the true Gamow-Teller transitions. 

It is worth mentioning that, recent studies of single 
and double beta-decays as well as of neutrino-nucleus re¬ 
actions under stellar conditions, have demonstrated an 
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TABLE III: The shift of the spectrum seperately of each state in MeV 


Positive Parity States Negative Parity States 


J+ 


32g 

48 




J- 


32^ 

48 Ti 




0 + 

2.60 

0.00 

0.65 

1.60 

0.90 

1.00 

0 " 

4.20 

1.00 

4.00 

4.30 

5.00 

4.47 

1 + 

5.00 

2.50 

2.65 

5.90 

2.50 

2.85 

1 ~ 

4.40 

4.05 

4.00 

4.20 

6.80 

4.30 

2 + 

4.35 

2.43 

2.10 

3.10 

2.55 

2.78 

2 " 

5.80 

4.40 

5.10 

6.80 

3.85 

2.39 

3+ 

5.90 

0.00 

2.70 

2.30 

2.50 

2.82 

3- 

6.00 

3.98 

4.10 

6.80 

2.60 

2.59 

4+ 

4.90 

3.56 

3.25 

2.50 

1.75 

0.00 

4~ 

5.00 

2.57 

4.25 

3.50 

3.55 

1.30 

5+ 

2.70 

0.84 

3.35 

2.00 

0.55 

2.40 

5" 

6.50 

0.00 

3.05 

3.50 

3.00 

0.00 


3.0 


2.0 


1.0 


0.0 


QRPA exp. 




FIG. 1: Comparison of the theoretical excitation spectrum, resulting from the solution of the pn-QRPA eigenvalue problem, 
with the low-lying (up to about 3 MeV) experimental one for and nuclei (for the other spectra see Ref. |29l I62 |L 

The agreement is quite good at least for low excitation energies. 


important role of the quenched value of the coupling con¬ 
stant gA El EH Eg. In the present calculations we also 
use a quenched value of gA same for all multipole transi¬ 
tions (see the following section). 


IV. RESULTS AND DISCUSSION 

In the ordinary muon capture on complex (A > 12) 
nuclei, the nuclear response is governed by the momen¬ 
tum transfer q of Eq. @, i.e. by an energy transfer to 
the daughter nucleus of the order of the muon mass 


minus the binding energy eb of the muon in the muonic 
atom restricted from below by the mass difference of the 
initial and final nuclei and from above by the muon mass 
[see Eq. @]. The phase space and the nuclear response 
favor lower nuclear excitations, namely the nuclear states 
in the giant resonance region (GDR and GT resonance) 
are expected to dominate [8j. 

In our calculational procedure we followed three steps, 
(i) In the first step we performed realistic state-by- 
state calculations on exclusive OMG rates in the isotopes 
Fe,®® Zn and ^°Zr, a set which covers 
a rather wide range of the periodic Table from light- to 





















































































































































7 


medium-weight nuclei. These calculations have been per¬ 
formed twice: Once with the use of the free nucleon cou¬ 
pling constants gA = 1.262 and the other with the use 
of the value gA = 1.135, to take into acount the rather 
small quenching effect indicated for medium-weight nu¬ 
clei [23l [53l |54] . We also focused on the study of the rela¬ 
tive strength of the polar-vector and axial-vector contri¬ 
butions for each individual excitation induced by the re¬ 
spective components of the muon-capture operators, (ii) 
In the second step of our calculations, we examined the 
dominance of the low-spin multipolarities into the total 
/r-capture rate. We also estimated the percentage (por¬ 
tion) of their contribution in the total rate for the most 
important multipolarities, (iii) In the last step, we eval¬ 
uated total muon-capture rates for the above set of iso¬ 
topes. For all the above calculations, the required wave 
functions (for the initial (ground) state and for all acces¬ 
sible final states) were constructed by solving the BCS 
and QRPA equati ons, respe ctively, as described before 
(see Sections III A and III BI. 


A. State by State calculations of exclusive 
transition Rates in p-capture 


At first, we evaluated the exclusive g -capture rates 
of Eq § for all multipolarities with J'^ < 5^. 
In Eq. (§ transitions between the ground state |i) = 
| 0 +) of a spherical target nucleus and an excited state 
IJJ) = I/) of the resulting odd-odd nucleus are consid¬ 
ered. In most of the previous studies a mean value of the 
muon wave function, with being the spheri¬ 

cal coordinate, has been utilized (see Appendix]^. An 
accurate description of the reaction Q (and of any re¬ 
action having the same initial state with it, i.e. a muon 
orbiting around an atomic nucleus {A,Z)), however, re¬ 
quires the exact muon wave function derived by solving 
the Schrodinger equation (or the Dirac equations) that 
obeys a bound muon within the extended Coulomb field 
of the nucleus in such muonic atoms |23| . 

Assuming, that the muon wave function in the region 
of the nuclear target is nearly constant, the integrals en¬ 
tering Eq. (§ can be performed by taking out of them an 
average value ($ 15 )■ Hence, the exclusive muon capture 
rates Ajj can be rewritten as: 


Ags^JJ - AjJ 


2G^{^uyRfgj- 

\{JJUMj - Cj)\\0j,) 


|(^/ll(r/-77“mio+)|^ (16) 


On the basis of the latter expression, we initially, per¬ 
formed state-by-state calculations, for the above men¬ 
tioned set of nuclear isotopes, by using the free nucleon 
coupling constant gA for the axial-vector form factor. 
Then, we repeated these calculations (with the excep¬ 
tion of and isotopes) by taking into account 
the quenching effect of the axial-vector coupling constant 


gA = 1.135. For each excitation of the daughter nucleus, 
our code provides us with the separate contributions in¬ 
duced by the components of the muon-capture operator. 
Relying on this possibility, we examined the multipole de¬ 
composition of the QRPA response in the muon capture 
reaction for the studied nuclei. In Figs. M and 1 ^ we 
illustrate the contribution of each individual transition. 
We also show the contribution of the polar-vector as well 
as the axial-vector parts originated from the correspond¬ 
ing components of the weak interaction Hamiltonian (see 
Sect. [h|. Evidently, most of the muon capture strength 
goes to I“, 1 ^ and 2 “ low-lying multipole excitations, 
of the particle bound spectrum and of the giant dipole, 
spin and spin-dipole resonances. 

As mentioned before, our code initially gives results 
for exclusive muon capture rates, Ajj, seperately for 
each multipolarity (in ascending order with respect to 
the pn-QRPA excitation energy Hj„). In order to study 
the dependence of the rates on the excitation energy lu 
throughout the entire pn-QRPA spectrum of the daugh¬ 
ter isotopes, a rearrangement of all possible excitations 
in ascending order with respect to w and with the cor¬ 
responding rates, is required. This was performed by 
using a special code (appropriate for matrices). Totally, 
for < 5^ in the model space chosen for each isotope, 
we have 286 states for the isotope, 440 states for 
each of the and '^^Ti isotopes, 488 states for each 
of the and and 912 states for isotope 

in the corresponding daughter nucleus. The variation of 
the exclusive rates throughout the entire excitation spec¬ 
trum of the daughter nucleus in the case of the above 
target isotopes are demonstrated in Figs. and 13 

For all reactions, the rates present some characteristic 
clearly pronounced peaks at various excitation energies 
uj and specifically for transitions J'^ = I'*', 1 “ but also 
for = O'*", 0“ and 2“ transitions. 

More specifically, in the daughter isotope the 

maximum peak corresponds to the I 7 QRPA transition 
at w = 7.712 MeV (see Fig. [ 3 . Other two characteristic 
peaks are at w = 18.135 MeH and at oj = 18.261 MeE 
which correspond to the Og and l^g transitions respec¬ 
tively. In the case of isotope the maximum peak 
corresponds to the 1^ transition at w = 4.855 MeE. An¬ 
other characteristic peak is at w = 15.564 MeE which 
corresponds to the Igg transition as shown in Fig. [^(left). 
For the '^^Sc isotope, the pronounced peaks correspond 
to the first excited 0+ state (0)^) (at w = 4.319MeE), 
the 2)"7 (at uj = 9.672 MeV), the l^^'g (w = 10.666MeE) 
and the l^g transitions (w = 18.868MeE). From Fig. 
13 (right pannel), for the daughter isotope ®®Mn, we see 
that the maximum peak appears at w = 8.278 MeE and 
corresponds to the I)*)) transition. Another important 
transition is that of Igg at uj = 18.716 A/eE. As shown 
in Fig 13 in the case of the daughter isotope ^^Cu, the 
maximum peak appears at w = 6.555 MeE and cor¬ 
responds to l)*g state and a pronounced peak for the 
l;^g at w = 14.833 MeE. Finally, for the isotope, 
the maximum peak appears for the 1^4 transition at 
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FIG. 2: Individual contribution of the Polar-Vector Av (pannel(a)) and Axial-Vector AA(pannel (b)) to the total muon-capture 
rate (pannel(c)) as a function of the excitation energy to for the and nuclei. 


w = 18.218 Mel^ and for the l^g at w = 9.752 MeF. 

From the above results, we conclude that in general, 
a great part of the OMC rate comes from the excitation 
energy region where the centroid of the GT strength is 
located for each daughter nucleus. As it is known from 
closure approximation studies [13 [H], the mean excita¬ 
tion energy in muon capture (about 15 MeV) is nearly 
equal to the energy of the giant dipole resonance (GDR) 
which is slowly decreasing with A or Z . On the other 
hand, the GT-like operators (in which the full spherical 
Bessel functions is taken into account) contribute very 
little in heavier nuclei where most of the active neutrons 
and protons are in different oscillator shells. In lighter nu¬ 
clei, however, i.e. for nuclei having N and Z smaller than 
40, the GT strength is significant and it is concentrated 
at the low energy region. Regarding the giant spin reso¬ 
nance (J'^ = I’*') for all nuclei the peak of the exclusive 
/i-capture rate is located between 5-11 MeV. It should 


be stressed that concerning the pronounced contribution 
to the 1“ states, it may contain a small portion of the 
spurious center of mass motion part (up to about 17% 
in our QRPA method) [^. This is due to the isoscalar 
movement of the nucleons in the mean field (dipole os¬ 
cillation of the whole nucleus). As it is known this is 
usually removed by using specific methods [55] . 

As it becomes clear from Figs. ii andj^ for the stud¬ 
ied nuclei the muon capture response presents maximum 
peak in the very important giant dipole resonance region 
(GDR), which is located in the energy region of 18-19 
MeV for and isotopes, and in the 

region of 15-16 MeV for and isotopes. These re¬ 
sults can be compared with the empirical expression, for 
medium-weight and heavy isotopes, which gives the en¬ 
ergy location of the giant dipole resonance, Eiyjj, based 
on the Jensen-Steinwedel and Goldhaber-Teller models 












































9 



. . ' .. '. ', , 

0 5 10 15 20 25 30 


a;(MeV) 



j rrT, ,^i ! i', iVr, , , , i , , , , 

0 5 10 15 20 25 30 


a;(MeV) 


FIG. 3: The same as Fig. [^but for the nuclei and ^^Fe. 


(a hydrodynamical view of the giant resonance) as |66| 

Eivd = 31.2A-1/3 + 20.6A-1/® (17) 

(A is the atomic mass of the nucleus). Even though this 
formula refers to pp- and nn-reactions, it can however be 
used to our results referred to pn-reactions (/r“-capture), 
on the basis of the well known Foldy-Walecka theorem 
according to which the giant dipole resonance in /re¬ 
capture rates, are calculated starting from the experi¬ 
mental phot o-absorption cross sections m- According 
to Eq. ( [TtI ) for the maximum 1“ peak is located 
at 18.668 MeV, for at 18.716 MeV, for at 

17.945 MeV and for at 16.684 MeV which are in 

a good agreement with our results (the worst case oc¬ 
curs for where the empirical peak is at about 15 

MeV). Moreover, our results are in good agreement with 
the conclusions of Ref. [24] where authors mention that 
for the stable Ni isotopes (58.60,62^^^ peak appears 
in the range of 18-19 MeV. We note that similar conclu¬ 


sion is extracted from the study of the charged current 
reaction ®®Fe(j^e, e“)^®C'o by Kolbe and Langanke [27] . 
where the peak of the giant dipole resonance appears at 
about 17 MeV (see Fig 1 of Ref. [17]). 

As can be seen from Figs. Hi and [^ the main con¬ 
tributions coming from the polar-vector operator are the 
1 “ and 0^ states while, the most important transitions 
due to the axial-vector operator are the 0“, 1+ and 2“ 
excitations, namely the lowest spin states. 

We note that the figures of this section, have been de¬ 
signed by using the ROOT program of Cern with binning 
width 0.112, 0.105, 0.105, 0.15, 0.14, 0.11, respectively, 
for ^^Zn, ^^Zr nuclei. 


B. Contribution of Multipole Transitions 

The second step of our study includes calculations of 
the partial /r“-capture rates for various low-spin multi- 
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FIG. 4: The same as Fig. [^but for the nuclei and 


polarities, Aj^r (for < 4=*=), in the chosen set of nuclei. 
These partial rates have been found by summing over the 
contibutions of all the individual multipole states of the 
studied multipolarity as 

Aj. = ( 18 ) 

/ 

/ 

+ -r7“"”)llo+)|"' 

/ 

where / runs over all states of the multipolarity |J^). 
As mentioned before, these calculations have been per¬ 
formed first by using the free nucleon axial-vector cou¬ 
pling constant qa = 1.262, and then by taking into ac¬ 
count the quenching effect indicated for medium-weight 
nuclei with gA = 1.135. 

For the target [STUSS] and isotopes, these cal¬ 


culations were performed only for the free nucleon cou¬ 
pling constant gA = 1.262. (the quenching effect can be 
ignored [25]). The results obtained for the partial /re¬ 
capture rates of these isotopes are illustrated in Fig. 
from which one can see that, as it is expected, the most 
important multipole transitions are the J'^ = 1'*" and 1~. 
More specifically, for isotope, the contributions of all 
= 1“ transitions exhaust the 36% of the total muon- 
capture rate and the = 1+ about 30%. Significant 
contribution, about 14%, comes from the J'^ = 0“ mul¬ 
tipolarity and about 13% from the J'^ = 2~. A similar 
picture is found in ^^5 isotope, where the dominant con¬ 
tributions to the total muon-capture rate are the = 1~ 
(38%) and the J'^ = 1+ (30%). From the rest of the 
multipolarities rather significant portions come from the 
abnormal parity transitions 0“ and 2“ about 13% and 
14% respectively. 

Because, as mentioned in the Indroduction, for elec¬ 
tromagnetic and weak charged-current nuclear processes. 
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FIG. 5: Partial muon capture rates Aj^r of different multipole transitions in and ^^5" isotopes. In both isotopes the 
pronounced contributions are the J’^ = 1“ and = 1^ multipolarity. 


the free nucleon coupling constant qa must be modified 
for medium-weight and heavy nuclei |23| . in /^“-capture 
on Fe,®® Zn,®® Zr isotopes we repeated the state- 

by-state calculations by using gA = 1.135 (a value smaller 
by about 10— 12% compared to the gA = 1.262). Histor¬ 
ically, the necessity of the renormalization of gA, came 
out of the following studies: (i) In the analysis of mea¬ 
surements on the nuclear beta-decays that lead to low- 
lying excitations |54j , and (ii) in the interpretation of the 
missing Gamow-Teller strength revealed in forward angle 
(p,n) and (n,p) charge-exchange reactions [53]. We note 
that, in (n,p) reactions many authors use quenched values 
of gA lying in the region of 0.9 < < 1-0 for nuclei with 

mass number 41 < H < 64 issiiMiizni- In /3 -decay and 
(p,n) reactions the quenching is mainly related to the ne¬ 
glect of configurations outside the model space used and 
the non-consideration of the meson-exchange currents. 


A quenched value of gA was recently suggested to be 
used in other weak interaction processes such as the neu¬ 
trino induced nuclear reactions. As has been found |25j . 
the consideration of a quenched factor instead of the free 
nucleon axial-vector coupling constant, leads to better 
agreement of the theoretical results with the experimen¬ 
tal muon capture rates. Since the axial-vector form fact or 
FA{q^) multiplies all four operators (see Eqs. (Al |-(A4|), 
a quenched value of gA must enter the multipole opera¬ 
tors generating the pronounced excitations 0“, 1^... etc. 
In Ref. [23] , a quenched value of gA is used only for the 
true Gamow-Teller transitions. In our study, we find that 
for the reproducibility of the experimental data, as the 
mass number A of the nucleus increases the quenching 
becomes more signifhcant and can not be ignored as we 
have done in the case of the ^*51 and isotope. 

For the medium-weight nuclei Zn and 

^°Zr, we used the moderate quenched value gA = 1.135 
and found that our rates are in good agreement with 
the results of other works [23]. By using this value of 
gA for the contributions of the different multipole transi¬ 
tions in the isotopes Zn and we found that 

the most important peaks correspond to the = I’*' 


and 1“. For the isotope, however, we found that a 
great part of the total rate comes from the = 1“ and 
2“, as is shown in Fig. In more detail, in the case 
of isotope the 1“ multipolarity contributes about 
44%, the 2“ about 17%, the 1+ about 16% and the 0“ 
about 11%. Significant contribution (about 7%) origi¬ 
nates also from the 0'*' multipolarity. For ®®Fe isotope 
the most important contribution about 42% comes from 
the 1“ multipolarity. Other multipolarities with signif¬ 
icant contributions are the l+(22%), 2“(13%), 0”(10%) 
and 0+(8%). A similar picture appears in the other two 
isotopes, Zn and Zr, where the major contribution is 
derived from the 1“ multipolarity, about 44% and 42% 
respectevelly. The 1“'" multipolarity contributes about 
21% in and about 20% in isotope. Corre¬ 

spondingly, the 2“ contributes about 13% for ®®Zn and 
about 14% for the 0+ about 8% and 9% respectiv- 
elly and finally the 0“ multipolarity offers about 8% for 
®®Zn and about 7% for 


In Table IV we present the partial muon-capture rates 
obtained for the low-spin multipole transitions up to 
J'^ = 4^ evaluated with our pn-QRPA code. Corre¬ 
spondingly, in Table [V] we tabulate the individual por¬ 
tions to the total OMC rate, for the low-spin nultipole 
transitions up to = 4^. As can be seen, for all nu¬ 
clei the contribution of 1“ multipole transitions is the 
most important multipolarity, exhausting more than 39% 
of the total muon-capture rate. Ordinary muon capture 
proceeds mainly through spin-multipole transitions, the 
most important of which are the Gamow-Teller transi¬ 
tions {j(i(kr)(jt'^ operator), and the spin-dipole transi¬ 
tions {ji{kr)\Yi ® operator) where jo and ji are 

the spherical Bessel functions of zero and first order, re¬ 
spectively |24j . Such important contribution is found in 
and in isotopes studied in Ref. HO]. 


There are no similar results for the isotopes 
^^Ti, ®®Fe and to compare with our portions. For 

the ^'^Zr, however, Kolbe, Langanke and Vogel [10] found 
about 28% (for 1“), 25% (for 1+) and about 13% (for 
2~) multipolarities which, with the exception of 1“ con- 
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FIG. 6: Contribution of multipole transition rates Aj^r (up to =4^) with the total muon capture rate in *^Ti, ^^Fe, 
Zn and ®°Zr isotopes with (filled histograms) and without (double dashed histogramms) quenching effect. The dominance 
of = 1“ and T*" multipolarities is obvious in all nuclei. 


TABLE IV: Muon capture rates Aj^r (in 10® s of each 
multipolarity evaluated with our pn-QRPA code. 




32^ 





0 " 

0.125 

0.168 

0.264 

0.398 

0.471 

0.662 

0 + 

0.037 

0.016 

0.177 

0.327 

0.488 

0.866 

1 " 

0.319 

0.481 

1.074 

1.740 

2.623 

4.087 

1 + 

0.271 

0.383 

0.397 

0.926 

1.263 

1.968 

2 " 

0.114 

0.171 

0.415 

0.524 

0.790 

1.307 

2 + 

0.014 

0.030 

0.065 

0.115 

0.169 

0.401 

3" 

0.001 

0.002 

0.006 

0.013 

0.020 

0.050 

3+ 

0.010 

0.012 

0.045 

0.073 

0.093 

0.255 

4" 

0.001 

0.001 

0.003 

0.006 

0.014 

0.029 

4+ 

0 .210"^ 

0.810"^ 

0 .210"® 

0.510"® 

0.710"® 

2.510"® 


TABLE V: The percentage of each multipolarity into the 
total muon-capture rate evaluated with our pn-QRPA code. 




32^ 

48yj 

sepe 



0 " 

14.03 

13.30 

10.78 

9.64 

7.94 

6.89 

0 + 

4.11 

1.27 

7.24 

7.92 

8.22 

8.99 

1 " 

35.74 

38.01 

43.88 

42.18 

44.21 

42.43 

1 + 

30.42 

30.28 

16.24 

22.46 

21.29 

20.43 

2 " 

12.81 

13.54 

16.97 

12.72 

13.32 

13.57 

2 + 

1.62 

2.36 

2.67 

2.79 

2.85 

4.16 

3" 

0.10 

0.15 

0.23 

0.32 

0.34 

0.52 

3+ 

1.09 

0.97 

1.82 

1.78 

1.58 

2.65 

4" 

0.06 

0.10 

0.14 

0.16 

0.23 

0.30 

4+ 

0.01 

0.01 

0.01 

0.01 

0.01 

0.03 


tribution, are in good agreement with our results listed 
in Table The difference in 1“ multipolarity is mostly 
due to the fact that is a double closed shell nucleus 
and the QRPA convergence is treated as in Ref. [Miin]- 


C. Total Muon-Capture-Rates 

In the last stage of our present work, we computed the 
total rates of muon-capture on the chosen set of nuclei. 
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These rates are obtained by summing over all partial mul¬ 
tipole transition rates in two steps. At first, we sum up 
the contribution of each final state of a specific multipo¬ 
larity, and then, we sum over the multipole responses (up 
to = 4^) as 

^tot = (19) 

j- j- / 

Such calculations have been carried out twice: one with 
gA = 1.262 (free nucleon axial-vector coupling constant) 
and the other with the quenched value gA = 1.135 [53]. 
The results are tabulated in Table [V^ where for the sake 
of comparisson we also include the experimental total 
rates as well as the theoretical ones of Ref. |53|. More¬ 
over, in Table [Vl| we show the individual contribution in 
the total muon capture rate of the polar-vector (A)(,j), 
the axial-vector (A(),(), and the overlap (A)^^) parts. 


calculated rates Zinner [53] obtained by using different 
approach and the comparison is good. 



TABLE VI: Individual contribution of Polar-vector, Axial- 
vector and Overlap part to the total muon-capture rate. Com¬ 
parison between the total muon captnre rates obtained by ns- 
ing the pn-QRPA with the quenched value of qa = 1.135 for 
medium-weight nucleus Ae,®® Zn and Zr) and the 

free nucleon coupling constant qa = 1.262 for the light nu¬ 
cleus Si and with the available experimental data and 
with the theoretical rates of Ref [23] . 



Total Muon-capture rates Atot(xl0®)s 

-1 


pn-QRPA Calculations 

Experiment RPA Dean 

Nucleus 

aV aa kVA a 

^Hot ^Hot ^Hot 

ALT 

ALT’- m 


0.150 0.751 -0.009 0.892 

0.871 

0.823 

32^ 

0.204 1.078 -0.017 1.265 

1.352 

1.269 

48 Ti 

0.628 1.902 -0.081 2.447 

2.590 

2.214 


1.075 3.179 -0.129 4.125 

4.411 

4.457 

^^Zn 

1.651 4.487 -0.204 5.934 

5.809 

4.976 

^°Zr 

2.679 7.310 -0.357 9.631 

9.350 

8.974 


As can be seen, our results obtained with the quenched 
gA are in very good agreement with the experimental to¬ 
tal muon-capture rates. For all studied nuclei the de¬ 
viations from the corresponding experimental rates are 
smaller than 7% when using the quenched gA (the devi¬ 
ation is much bigger when using the gA = 1.262). So, for 
the reliability of our results it is necessary to take into ac¬ 
count the quenching effect. To make it more perceptual, 
in Fig. [^we have plotted the ratio of our theoretical total 
muon-capture rates divited by the experimental ones, i.e. 


for the results obtained with the above two values of gA 
(with and without quenching). The filled circles repre¬ 
sent the results for the free gA and the X symbols the 
results for the quenched gA- It is evident the better 
agreement of our calculations with quenched value of gA- 
We furthermore, compare our results with the available 


FIG. 7: Ratio of the calculated and experimental total muon 
capture rates as a function of Z. Circles and X symbols corre¬ 
spond to rates calculated with free nucleon qa and quenched 
value of QA respectively. 

Finally, it is worth noticing that, in medium-weight 
nuclei the contribution comes mainly from transitions 
for which the angular momentum transfer is L=0,1 and 
2 but, in heavy nuclei, some contributions from higher 
multipolarities become noticeable. 

V. SUMMARY AND CONCLUSION 

In the present work, relying on an advantageous nu¬ 
merical approach constructed by our group recently, we 
performed detailed calculations for all multipole transi¬ 
tion matrix elements entering the exclusive muon-capture 
rates. The required nuclear wave functions were obtained 
within the context of the pn-QRPA using realistic two- 
body forces (Bonn C-D potential). Results for the exclu¬ 
sive rates through extensive state-by-state calculations 
and subsequently for the total muon capture rates on 
the set of isotopes and 

were computed. 

Because the capture rates are rather sensitive to the 
quenching of the axial-vector coupling constant, we ex¬ 
amined the known as in-medium effect of the nucleon, 
by reducing this constant from its free nucleon value 
gA = 1.262 to the effective value gA = 1.135 for all multi¬ 
pole transitions, and found that the experimental muon 
capture rates are well reproduced with an accuracy bet¬ 
ter than 10%. Detailed study of this effect, however, 
required for experiments at RCNP [71175] is under way 
and results are expected to be obtained soon. 

The muon-capture studies on these nuclei demonstrate 
that the used pn-QRPA method may provide an accurate 
description of the charged current semileptonic weak in¬ 
teraction processes in the Z-range of the isotopes chosen. 
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As the inclusive muon capture rates and the cross section 
of the antineutrino-induced charged current reactions are 
closely related (both of them are governed by the same 
nuclear matrix elements and proceed via the same initial 
and final states), we have adopted this method to study 
other types of charge-changing weak interaction processes 
as, electron-capture, beta-decay modes, etc. ESIISI] in 
currently interesting nuclei from a nuclear astrophysics 
point of view. 
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Appendix A: Nuclear Matrix Elements 

The ei gh^ t different tensor multipole operators enter¬ 
ing Eq. ^ refer to as Coulomb Mj, longitudinal Cj, 
transverse electric Tf and transverse magnetic 
contain polar-vector as well as axial-vector parts and are 
written as: 


elements of the form (ji||T/||j2), where T/ represents 
any of the seven basic multipole oper a tors {M-h, 


^jv/’ ^m) ( 

AI 

)-( 

A4 

), have been 

written in closed compact formulae 

as [341 li 

M 


n-max 

{{n,h)n\\T-^\\in2l2)j2) = e-yyP/^ £ (A5) 

fi—0 


where the coefficients are given in Ref. [3l]. In the 
latter summation the upper index n^ax represents the 
maximun harmonic oscillator quanta included in the ac¬ 
tive model space chosen as Umax = (A^i + N 2 — /3)/2, 
where Ni = 2ni + /i, i=I,2, and /3 is related to the rank 
of the above operators [33]. 

In the context of the pn-QRPA, the required reduced 
nuclear matrix element between the initial |0g^) and any 
final I/) state entering the rates of Eq. ( [l^ are given by 

(/iif^iioj) = i: 

j2>jl 

■ (A6) 

where Uj and Vj are the probability amplitudes for the j- 
level to be unoccupied or occupied, respectively (see the 
text) m- 

These matrix elements enter the description of various 
semi-leptonic weak interaction processes in the presence 

of nuclei [31111USHISIIMHIS]. 


MjM{qr) = (AI) 

= Pi - i^[FAni^{qr) 

+ 2^P^ P qoFp)'P‘M ( 9 ^)] 

PjAtiq^) = Ljm P L^jm (A2) 

= '^FYMij{qr)+iFAT:P{qr) 

PjM^q"^) = PjM + PjM (A3) 

= ^Miqr) + ^^J-^yM{qr)] 

+ iFAFj'Yiqr) 

_ i^magn r^magnb 

- \f^^y'M{qi")]+iFAyM{qr) (A4) 


where the form factors Fx, X=I,A,P and /i'^ are func¬ 
tions of the 4-momentum transfer . 

J-fl 

These multipole operators, due to the Conserved Vec¬ 
tor Current (CVC) theory, are reduced to seven new ba¬ 
sic operators expressed in terms of spherical Bessel func¬ 
tions, spherical harmonics and vector spherical harmon¬ 
ics (see Refs. |il[53]). The single particle reduced matrix 


Appendix B: Muon wave function in the muonic 
atom 

The calculation of the exact muon wave function, 
$is(r), entering Eq. ^ needs the use of a specific nu¬ 
merical method. This however, can be avoided by using 
either its value at r ~ 0, namely the 4’is(r ~ 0), or as 
stated in Sect. |llj an average value (‘I’ls), which is given 
in terms of the effective nuclear charge Zef / that sees the 
bound muon as 

= (Bl) 

TT ^ Z 

{a denotes the fine structure constant). The quantity 
Zgf f is approximated by Z^^ = ttcKq (p), where ag is the 
muon Bohr radius and (p) is the mean charge density 
of the parent nucleus m- For light nuclei Z^ff — Z 
but for heavier ones Z^ff ^ Z. In recent studies the 
exact wave functions for the bound muon are obtained 
by solving the Schroedinger and Dirac equations by using 
neural network techniques or genetic algorithms (73). In 
the work of Zinner, Langanke and Vogel [23], for the 
description of the exact bound muon wave functions (w- 
fs), the muon density beyond the site of the nucleus is 
considered for solving the Dirac equation. These authors 
use exact muon wave functions, for other muonic orbits, 
$2p, etc, which are considered to have rather signifhcant 
contributions [23] . 
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Appendix C: pn-QRPA equations 


The new amplitudes satisfy the matrix expressions 


In our numerical solution performance we rewrite the 
QRPA equations (15) by defining a new set of ampli¬ 
tudes P™ and i?"* which are related to the preceding 
ones through 


y™ = (Cl) 


(A - B)P^ = P™ , (A -f P)P™ = (C2) 


Then, we have 

iA + B){A-B)P^ = nl^P"^. (C3) 

The latter equations can be diagonalized seperately and, 
subsequently, the X and Y amplitudes are directly deter¬ 
mined [ISJI3S] . 
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